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We discuss heat conductivity from the point of view of a variational multi-fluid 
model, treating entropy as a dynamical entity. We demonstrate that a two-fluid 
model with a massive fluid component and a massless entropy can reproduce a 
number of key results from extended irreversible thermodynamics. In particular, we 
show that the entropy entrainment is intimately linked to the thermal relaxation 
time that is required to make heat propagation in solids causal. We also discuss 
non-local terms that arise naturally in a dissipative multi-fluid model, and relate 
these terms to those of phonon hydrodynamics. Finally, we formulate a complete 
heat conducting two-component model and discuss briefly the new dissipative terms 
that arise. 
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1. Introduction 

Heat conductivity is a central problem in thermodynamics. It is well known that 
the classical description in irreversible thermodynamics, essentially Fourier's law, 
has unattractive features. In particular, it predicts an instantaneous propagation of 
thermal signals. This is in contradiction to the expected hyperbolic nature of phys- 
ical laws. In fact, the associated non-causality would be completely unacceptable 
within a relativistic theory. Resolving this issue has been a main motivating factor 
behind the development of extended irreversible thermodynamics (Jou et al, 1993; 
Miiller & Ruggeri, 1993), a model which introduces additional dynamical flelds in 
order to retain hyperbolicity and causality. This model has proved useful in various 
application areas, ranging from superfluid systems to heat conduction in solids. 

In this paper we consider the problem of heat conductivity within the flux- 
conservative multi-fluid framework developed by Andersson & Comer (2006) (with 
the corrections discussed by Haskell et al (2009)). This model builds on the non- 
dissipative variational model developed by Prix (2004), and represents the natu- 
ral non-relativistic counterpart to Carter's convective variational hydrodynamics in 
general relativity (Carter, 1989), see Andersson & Comer (2007) for a recent review. 
These models have so far primarily been used to investigate the dynamics of com- 
pact stars (see Andersson & Comer (2001); Glampedakis et al (2007); Glampedakis 
& Andersson (2009); Passamonti et al (2009) for details). Yet, one would expect 
the general framework to be universally relevant. It is therefore of some interest 
to consider applications in other problem areas. This serves several purposes. First 
of all, it is often the case that the elegant geometric view of a relativistic analysis 
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(where time and space are treated "equally" ) simplifies the description of a complex 
system. The close link between the variational multi-fluid model and the relativistic 
counterpart may lead to insights that would be hard to reach otherwise. Secondly, 
it is interesting to learn from the requirements of models for a range of different 
systems. Ultimately, this will improve our understanding of the general multi- fluid 
framework, and the role of the various parameters in the model. 

In this paper we focus on the simplest "conducting" system, with a single species 
of particle together with a massless entropy component. If the entropy drifts relative 
to the particles, the system is heat conducting. Throughout the discussion we will 
assume that the entropy can be treated as a (massless) fluid. In essence, this means 
that the "phonon mean free path" is not too large. We aim to demonstrate that 
our flux-conservative formulation for this system captures key aspects of extended 
thermodynamics. We illustrate this by writing down a model for heat conductivity 
in a rigid solid and comparing to various results in the literature. This exercise 
makes it clear that our model incorporates a finite propagation speed for heat. We 
also learn that the associated relaxation timescale is directly linked to the entrain- 
ment between particles and entropy. This entropy entrainment also played a central 
role in our recent discussion of finite temperature supcrfiuids (Andcrsson & Comer, 
2009). In fact, the model we consider here is formally equivalent to our recent model 
for superfluid Helium, the key element being that the massless entropy component 
is allowed to flow relative to the particles in the system. However, the discussion 
here differs in that we do not impose the irrotationality constraint associated with 
a superfluid system. Instead, we focus on the heat conductivity. We wish to under- 
stand to what extent this, conceptually rather elegant, model captures the complex 
physics associated with heat flow. As an interesting by-product, the present discus- 
sion suggests how the Helium model could be extended (using a three-fluid model) 
to account for the thermal conductivity associated with the interaction between 
rotons and phonons that dominates at higher temperatures (Khalatnikov, 1965). 
This is an interesting problem since the development of a causal model for this 
phenomenon is still outstanding. Our analysis also hints at the key ingredients of a 
causal relativistic model for heat conductivity. The main lesson is that the entropy 
entrainment must be retained in order to avoid pathological behaviour. 

2. Multi-Fluid formulation 

We consider a simple system with two dynamical degrees of freedom. We distinguish 
between the mass carrying "atoms" and the massless "entropy" . The former will be 
identifled by constituent index n, it has mass m, number density n and flows with a 
velocity ?;* , while the latter is represented by s, with number density s (the entropy 
per unit volume) and a flow given by vl- The canonical momentum (density) for 
each fluid is determined by the variational analysis of Prix (2004). This leads to 

< = mnvf - 2aw';' , (2.1) 

for the atoms and 

TTf = 2awr , (2.2) 

for the entropy. This latter relation encodes the inertia of heat. We see that the 
heat's inertia vanishes if the entrainment between the two components, quantified 
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by a, vanishes. In these expressions we have used the relative velocity wf^ = — 
v^. We recall that the variational analysis takes as its starting point an energy 
functional E, representing the equation of state. The coefficient a then follows 
from 

dE 

a 



(2.3) 



As discussed by Andersson & Comer (2006, 2009) the momentum equations for 
the dissipative two-component system can be written 



and 



1 



(2.4) 



(2.5) 



Here the terms and D^j represent the dissipation and /* and /* are the respec- 
tive forces acting on the matter and entropy. The chemical potentials are determined 
by 

dE 



dn 



and 



Ms 



dE 
ds 



(2.6) 



(2.7) 



where we have identified the temperature, T, in the usual way (note that we use 
units such that Boltzmann's constant is equal to unity, ks = l)t- 

Since we have no particle creation or destruction, mass conservation leads to 



dtn + Vjinvi) =0 . 
At the same time the entropy can increase, so we have 



(2.8) 



(2.9) 



where the second law of thermodynamics requires Fg > 0. Finally, the dissipative 
terms are constrained by the fact that we consider a closed system. As discussed 
by Andersson & Comer (2006), this means that we must have 



f! 



-ft 



(2.10) 



and 

= A, - . (2.11) 

Here, Dij is the total "dissipation". It is distinguished by the fact that it is sym- 
metric in its indices whereas Df- and do not have to be. 

t It should be pointed out that the notion of temperature in a non-equiUbrium system is 
non-trivial, see Casas- Vazquez & Jou (2003) for a thorough discussion. In the present analysis we 
assume that the temperature is obtained from (2.7), i.e. in the same way as in thermal equilibrium. 
This "operational" definition seems the most natural in the present context. 
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Following the steps taken in Andersson & Comer (2009), i.e. constructing the 
dissipative terms from the relevant thermodynamic fluxes and taking account of 
the Onsager symmetry, we arrive at^ 

-/° = 27e""wr + 2<S"°Wf " , (2.12) 

-DJ'^. =5"%-,M;^, + 5,,(C"ens + Ce,)+277'^"e- + 2,y"e^,. + a'^"ei,fel^^ , (2.13) 
and 

-A, = (C"e„, + ce,) + 2?7'^e^; + 2r,e^^. . (2.14) 

Here we have defined the expansion 

Os = ^jvi , (2.15) 

the trace-free shear 

= I (Vi^,' + - IgijQ,^ , (2.16) 

and the "vorticity" 

W: = ^e^^''{V,vl-Vkv'^) , (2.17) 

of the entropy flow. The quantities 6ns, ©"/ and W^^ are constructed from w^^^ 
analogously. 

So far, the development has been formal. The model was developed by com- 
bining the non-dissipative equations of motion from the variational analysis (Prix, 
2004) with the general form for the dissipative terms assuming a quadratic deviation 
from thermodynamic equilibrium. We have accounted for the Onsager symmetry 
between the dissipative terms, and imposed the conservation of total angular mo- 
mentum. This reduces the problem to one with nine, as yet unspecified, dissipation 
coefficients; 71™, <S"", C", C™, C and ry. 

It should be noted that a number of, essentially, analogous formulations for 
multi-component systems exist in the literature, see for example Drew & Passman 
(2007); Rajagopal & Tao (1995). However, most of these arc not arrived at by varia- 
tional arguments, and hence they do not distinguish the fluxes from the canonically 
conjugate momenta. This means that the entrainment effect tends to be ignored. 
In fact, the novelty of our approach relates to the entrainment between particles 
and entropy. From an intuitive point of view, this encodes the inertia of heat and 
allows us to assign an effective mass to the entropy component. It is well known 
that these concepts are central to causal models of heat conductivity, and we will 
soon see why this is so. 

X We use a coordinate basis to represent tensorial relations. In other words, we distinguish 
between co- and contra- variant objects, Vi and v^, respectively. Indices, which range from 1 to 
3, can be raised and lowered with the (flat space) metric gij, i.e., vt = gijV^ . Derivatives are 
expressed in terms of the covariant derivative Vj which is consistent with the metric in the sense 
that Vjgfcj = 0. This formulation has great advantage when one wants to discuss the geometric 
nature of the different dissipation coefficients. We also use the volume form eijk which is completely 
antisymmetric, and which has only one independent component (equal to ^/p in the present case). 
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3. The standard approach to heat conductivity 

It is natural to begin the discussion of heat conductivity by sketching the classic 
approach to the problem. This provides a useful contrast to the more general model 
that we will develop. 

The following argument is more or less identical to that provided by Prix (2004). 
First we rewrite (2.9) as 

dts + Vi[svl, + swlJ=T, . (3.1) 

In general, Fg > if the system is out of equilibrium. Next we introduce the heat 
flux vector 

q' = sT< , (3.2) 



such that 



dts + V, 



(3.3) 



Keeping only the linear friction term associated with TZ™^, and using eq. (2.9) from 
Andersson & Comer (2009), i.e. 

TFs = -i-<, - D\Vivl - D^^Vjvt , (3.4) 

we arrive at 

F, = -^/>L = |7^-^ns = 27^"" (^)' > . (3.5) 

The second law of thermodynamics thus requires that 7^""^ > 0. In the following 
wc will often omit the subscripts on the resistivity coefficient 7?."". This makes 
the equations somewhat clearer and we do not believe that it should cause any 
confusion. 

Finally, if we consider the case of vanishing entrainment (cf. Prix (2004)) then 
it follows from (2.5) that 

ft = sWiT = -/f = 2nwr = -2^qi , (3.6) 

or 

1 s^T 

* = -2^^^^- ^'-'^ 
Comparing to Fourier's law we identify the thermal conductivity as 



2 n 



(3.8) 



This completes the traditional description of the heat conductivity problem. How- 
ever, for reasons that we have already discussed, this analysis is not entirely satisfac- 
tory. Fortunately, the multi-fluid approach advocated here allows for a more general 
view. This alternative description has a richer phenomenology and recovers a num- 
ber of features of extended irreversible thermodynamics (Jou et al, 1993; Miiller & 
Ruggeri, 1993). Illustrating this is one of the main purposes of our discussion. 
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4. The extended thermodynamics view 

The analysis in the previous section focussed on the entropy conservation law. The 
spirit of the discussion was very much that of classic irreversible thermodynamics. 
We will now depart from this view, and ask what we can learn by taking the two- 
fluid model, e.g. the entropy inertia, at face value. The massless entropy flow and 
the associated entrainment then play central roles. 

Let us consider the entropy momentum equation (2.5) a bit closer. We can 
rewrite this equation as 



1 . „\ . f I 



-TT\ 



+ ,sV,T + nl{y,vi) + V,< + V,D% = f! . (4.1) 

For simplicity, we flrst include only the resistive contribution to the dissipation, i.e., 
we set 

Dl^=Q, and = 27^<^ (4.2) 

as before. Then we have 

+ sViT + nU^jvi) + n^^Vivi = ^< , (4.3) 

where we recall the deflnition of the entropy momentum, (2.2). 

It turns out to be instructive to view this equation in the context of extended 
irreversible thermodynamics (Jou et al, 1993; Miiller & Ruggeri, 1993). To facilitate 
the comparison, let us focus on the simplified case of heat conduction in a rigid solid. 
In that case we have vl^ = and our momentum equation simplifies to 

dt^i - Vi - '^.^ V. (^TT^) + sV,T = ^< , (4.4) 

Moreover, as before it makes sense to introduce the heat fiux as 

wl^ = -^q' = . (4.5) 

™ sT 2a ' ^ ^ 

In the following we will often choose to keep the relative velocity as the main variable 
(we also omit the constituent indices, using = wl^ for clarity). This simplifies 
many of the equations as it suppresses factors of sT. Expressing all results in terms 
of the heat flux vector is trivial given (4.5). 
It follows that 

-dt {2awi) - Wj {2aw'wi) - aViV? + sVjT + 2Tlwi = . (4.6) 

This is one of our main results. 

Consider the limit of vanishing entrainment. Letting a ^ we are left with 

-. = -4v.T, (4.7) 
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or ^ 

ft = -^V,T=-KV,r, (4.8) 

where we have identified the thermal conductivity, k, as in the previous section. In 
the hmit of vanishing entrainment we recover Fourier's law, as expected. 

If wc instead linearise the equation (with respect to thermal equilibrium, = 0) 
we find that 

'di(^^q,^+sV,T+^qi = . (4.9) 

For constant parameters!, this can be written 

Tdtqi + qi = -KViT , (4.10) 
where we have introduced the thermal relaxation time, r, according to 



a 

'n 



(4.11) 



Eq. (4.10) is known as the Cattaneo equation (Cattaneo, 1948). It resolves the 
"paradox" associated with instantaneous propagation of heat predicted in classic 
irreversible thermodynamics, see Jou & Casas- Vazquez (1988) for a brief discus- 
sion. It also leads to the presence of a second sound in solids, an efi'ect that has 
been observed in laboratory experiments on dielectric crystals, see for example the 
discussion by Ruggeri et al (1996). As described by Jou et al (1993), the Catta- 
neo equation provided key stimulus for the development of extended irreversible 
thermodynamics. 

The above exercise shows that the variational multi-fluid formalism contains a 
key element of extended thermodynamics. Moreover, we have a natural interpreta- 
tion of the entropy entrainment, a. In terms of the relaxation time r we have 

That is, wc should expect to have a. < 0. 

Is it surprising that the Cattaneo equation can be deduced from the variational 
two-fluid model? Not really. In retrospect the result is, more or less, obvious. This 
is easy to sec if wc compare the nature of the variational energy fimctional E (Prix, 
2004) to the generalised entropy used in extended irreversible thermodynamics (Jou 
& Casas- Vazquez, 1988). In the variational case, we have E = E{n,s,w'^) and in 
the case of a solid (taking n to be constant) it follows that 

dE dE / \ ^ 

dE = —ds + TT^dw'^ = Tds + adw^ = Tds + ad( — ] . (4.13) 

OS ow"^ \sT/ 

Meanwhile, one of the tenets of extended irreversible thermodynamics is the gen- 
eralised entropy 

s{u,q) = So{u)-^mQ\ (4-14) 

t There is no physical reason why a/sT should be "constant". We only maice the assumption 
in order to facilitate a direct compajrison with the Cattaneo equation. The discussion of Morro 
& Ruggeri (1988) makes it quite clear that the parameter needs to be temperature dependent in 
order for this kind of model to be able to reproduce data from second sound experiments. 
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where u is the internal energy. As discussed by Alvarez et al (2008), this assumption 
is common to a number of alternative approaches to irreversible thermodynamics. 
The entropy satisfies the generalised Gibbs identity 

ds = ^du - Pdq^ . (4.15) 

It is now apparent that, once we identify u = E, (4.13) and (4.15) contain similar 
information. The key point is that the cncrgy/cntropy of the system depends on 
the heat flux. This is natural since the heat flux corresponds to the flow of energy 
relative to the matter current. The variational approach provides a slightly different 
perspective on the problem, but the predicted dynamics should be equivalent. 

Before moving on, it is worth discussing the role of the non-linear terms that 
were discarded in (4.9). The variational analysis naturally leads to the presence of 
quadratic terms in the heat flux. Due to their origin, i.e. the entropy momentum 
equation, these terms (essentially) take the same form as the non-linear terms in 
the standard Euler equation. In other words, they do not represent the most general 
non- linearities that one might envisage (in the dissipativc problem). In the interest 
of clarity, we have chosen not to compare the non-linear terms in our model to var- 
ious attempts at constructing non-linear heat-conducting models, see for example 
Morro & Ruggeri (1987); Ruggcri ct al (1996); Jou et al (2004); Lcbon ct al (2008b) 
and Llebot et al (1983). Such a comparison would obviously be interesting given 
that non-linearities are relevant for the development of both shocks and turbulence. 
However, out main initial aim is to establish the viability of the multi-fluids ap- 
proach to the heat problem. For this purpose, the evidence provided by the linear 
comparison should be adequate. 

5. Including "non-local" terms 

The model discussed in the previous section is the simplest in a hierarchy of pos- 
sible heat conductivity models. By relaxing the assumptions that led to (4.10) we 
can easily obtain more complicated models. Such models are interesting, because 
(4.10) may not provide a faithful representation of all relevant phenomena. It is 
obviously important to compare a more general model to analogous efforts in the 
literature. This comparison provides further insight into the usefulness of the multi- 
fluid framework and the interpretation of the various coefficients. 

So far we have neglected the viscous stresses in the entropy momentum equation. 
Let us now relax this assumption, i.e., account also for D^j. In order to simplify the 
analysis somewhat, we will still assume that .S''" = a'^"^ = 0. That is, we ignore the 
coupling to vorticity. 

We now need 

(5.1) 

Rewriting this expression in terms of the particle velocity v^^ and the relative fiow 
we have 

D!j = -9ij (C - Oe„ - 2{rj - r/")e^,. - ff,,Ce„, - 2^6^/ , (5.2) 
where we have introduced 

^ = 2C" - C"' - C, and f] = 2rf - rf"" - ri . (5.3) 
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Again focussing on the case of a rigid solid, the first two terms in (5.2) vanish 
and we also have 

Ons = -^iw' , (5.4) 

and 

If 9 .1 

(5.5) 



e"? = - 

'3 2 



2 , 



Adding the relevant terms to (4.6), taking ^ and r] to be constant, for simplicity, 
we have 

- dt{2awi) - Vj {2aw^Wi) - aViw'^ + sViT 

72„ - 1 



+ 27^^^;i + fjV^Wi + [C + -^vj ViV^w^ = . (5.6) 

Wc can simplify things further by i) assuming that all parameters arc constant, 
and ii) linearising in (the caveats regarding these assumptions remain as before) . 
This leads to the equation 

-^d,q, + sV^T + ^g, + A v^g, + ^ (C + \n) V,V,g^ = , (5.7) 

or, making use of the thermal relaxation time and the thermal conductivity, 

Tdtq^ + kV,T + ft + ^2 v^ft + (^l + 1^2^ V,V,g^ = . (5.8) 



We have also introduced the two lengthscales 

27^ 



(5-9) 



(5.10) 



and 

|2-X 

These can be taken to represent the mean-free path associated with the dissipative 
terms. 

This result can be directly compared to the "phonon hydrodynamics" model de- 
veloped by Guyer & Krumhansl (1966) (see Llebot et al (1983), and Cimmelli (2007) 
for alternative descriptions). Their model is the most celebrated attempt to account 
for non-local heat conduction effects. It accounts for interaction of phonons with 
each other and the lattice. Resistive terms are represented by r while momentum 
conserving interactions are associated with and Z^. Our model, (5.8), completely 
reproduces the Guyer & Krumhansl (1966) result, once we set 



ll = -ll . (5.11) 



This leads to 

TdtQi + KViT + qi + ll {V^qi + 2ViV,g^) = . (5.12) 

The usefulness of this result is due to the fact that it can be used both in the 
collision dominated and the ballistic phonon regime. In the former, the resistivity 
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dominates, the nonlocal terms can be neglected and heat propagates as waves. In 
the opposite regime, the momentum conserving interactions are dominant and we 
can neglect the thermal relaxation. In this regime, heat propagates by diffusion. 
The transition between these two extremes has recently been discussed by Vasquez 
& Markus (2009). 

Interestingly, the non-local heat conduction model may be useful in the descrip- 
tion of nano-size systems. If a system has characteristic size L, and /^/L ^ 1, then 
one would not necessarily expect a fluid model to apply. Nevertheless, Alvarez et 
al (2009) have argued that the expected behaviour of the thermal conductivity as 
the size of the system decreases (as discussed by Alvarez & Jou (2007) one would 
expect the "effective" conductivity to scale as L/lr,) can be reproduced from (5.8) 
provided that the appropriate slip condition for is applied at the boundaries. A 
key part of this analysis is the close analogy between (5.8) and the Navier-Stokes 
equation. In the latter case it is well-known that an applied non-slip condition at 
a surface leads to the formation of a viscous boundary layer that dominates the 
dissipation of the bulk flow. It appears that the heat problem is quite similar in the 
ballistic phonon regime, although the required slip condition is different. This is an 
interesting problem that requires more detailed study. 



At this point we have demonstrated that the multi-fluid formalism, with one fluid 
representing the massless entropy, reproduces a number of non-trivial results for 
heat conductivity in rigid solids. However, this is a simplified problem since one of 
the degrees of freedom in the system was "frozen" . In order to complete the model, 
we will now relax this assumption and allow ^ 0. This leads to a system of 
equations with interesting applications. In particular, one could imagine modelling 
systems with spatial transitions to superfluidity (as in a neutron star core). In one 
regime, the system would be dominated by resistivity, while the thermal relaxation 
timescale determines the dynamics elsewhere. A unified model for this problem 
could prove very useful indeed. 

Let us return to the general two-fluid model and focus on the matter degree 
of freedom. As in the case of superfluid Helium (Andersson & Comer, 2009), it is 
natural to work with the total momentum equation. By combining (2.4) and (2.5) 
we have 



6. The general two-component model 




(6.1) 



We can rewrite this using 



(6.2) 




(6.3) 



and 



,ns 



(6.4) 
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We also use the fact that the pressure foUows from (Andersson & Comer, 2009) 

ViP = nViiin + sViT - aViwl^ . (6.5) 
Combining these results we arrive at 

dtiK) + V; ipvlvf) - Vi{2awlwr) + V^p = -VjD\ . (6.6) 

Using our definition for the heat flux, (4.5), we see that if we linearise in Wi (or equiv- 

alently, qi) then the problem simplifies considerably. If we also use the continuity 
equation (2.8), we arrive at an equation that resembles the standard Navier-Stokes 
result (Landau & Lifshitz, 1959); 



p{dt + vlyi)vf + Vip-- 



The right-hand side is, however, different. Keeping 
find that 



(6.7) 

(as before) we 



V,(VX) 



where we have defined 



(6.6 



C = - C, 



and 



= r? -ri . 



(6.9) 



For simplicity, we have also assumed that the equilibrium configuration is uniform. 
This leads to the final result 



1 



1 



C + ifl]^i{Vjq')-^V'qi . (6.10) 



The first two terms on the right-hand side are familiar from the Navier-Stokes equa- 
tion (Landau & Lifshitz, 1959), but the last two terms are new. They represent the 
dissipative coupling between the total momentum and the heat flux. An interest- 
ing question concerns whether there are situations where these terms have decisive 
impact on the dynamics. Are there, for example, situations where the last term is 
similar in magnitude to the second term? 

To complete the model, we need the momentum equation for the heat flux. 
Starting from (4.1), linearising in the heat flux and using the results from the 
previous section, we flnd that 



{dt + vlVj)qi + -qi + qiVjvi + qjWivi 

T 



K„ nn sT 

--ViT 



1 

+ 2^ 



C+^fj]Vi{Vjq^)+vV^qi 



(6.11) 



Comparing to the corresponding equation for a rigid solid, we recognize several 
terms. Some dissipative terms are, however, new. 



Article submitted to Royal Society 



12 



N. Andersson} and G. L. Corner^ 



Before we conclude our discussion, it is worth making the foUowing observation. 
In our formulation of the problem we have assumed that the entropy conservation 
law is used explicitly. An alternative, more common, strategy is to use the energy 
equation. The two approaches are, in principle, equivalent. Nevertheless, it is useful 
to complement our analysis with a brief consideration of the energy equation. Prom 
Andersson & Comer (2006) we have (for an isolated system) the energy equation 

dtU + ViQ' = . (6.12) 

After some work, using the various definitions, we find that 

Q' = {^^Pvl - 2a< + n/x + sT^ < - {sT + 2a^;>7) . (6.13) 

Here, we can use the fundamental relation 

p + E = nn + sT, (6.14) 

to get 

Q'=[\pvl- 2a< +p + E^vi-{sT + 2at;>f ) wi, . (6.15) 
We also have 

U=^pvl-2awl + E . (6.16) 

Combining these results, and linearising in w^g, we immediately arrive at the stan- 
dard energy equation (c.f. Landau & Lifshitz (1959)) 

= Vi (sT<J = -Viq' . (6.17) 

This confirms, at least at the linear level, the definition (3.2) of the heat flux. 

We now have a "complete" model for a heat conducting fluid. It combines the 
equations of motion (6.10) and (6.11) with the two continuity equations (2.8) and 
(2.9). This model should be relevant for dynamics on timescales such that the ther- 
mal relaxation cannot be ignored. In a typical system, this would correspond to 
the extreme high-frequency regime. However, as wc already know from the discus- 
sion of Andersson & Comer (2009) the model also applies to superfluid condensates 
at finite temperatures. In essence, we have a unified framework for modelling the 
transition to superfluidity. Finally, there may be situations where the model ap- 
plies even if the thermal relaxation can be safely ignored. This would be the case 
when the "phonon" mean free path exceeds the size of the system. As discussed 
in the previous section, the dissipativc terms in (6.11) then play the leading role 
and the crucial boundary effects may be incorporated by imposing suitable surface 
conditions on the heat flux. This possibility has so far been discussed only for nano- 
systems, but it is worth noting that it may be relevant also for large scale systems. 
In fact, the ballistic phonon regime should apply to cold superfluid condensates in 
neutron stars. To what extent the present analysis can be applied to that problem 
is an interesting question for the future. 
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7. Concluding remarks 

We have discussed heat conductivity from the point of view of the variational multi- 
fluid model developed by Prix (2004) and Andersson & Comer (2006). We have 
shown that a two-fluid model that distinguishes between a massive fluid compo- 
nent and a massless entropy can reproduce a number of key results from extended 
irreversible thermodynamics. In particular, we have demonstrated that the entropy 
cntrainmcnt, that played a central role in our recent discussion of supcrfluicl He- 
lium (Andersson & Comer, 2009), is intimately linked to the thermal relaxation 
time that is required to make heat propagation in solids causal. We have also con- 
sidered non-local terms that arise naturally in the dissipative multi-fluid model, 
and related them to models of phonon hydrodynamics. This discussion may pro- 
vide useful insight into the modelling of both nano-systems and superfluids at low 
temperatures where the phonon mean-free path is large compared to the size of 
the system. Finally, we formulated a "complete" heat conducting two-fluid model 
and identified a number of "new" dissipative terms. Future work needs to establish 
whether there are physical situations where these terms play a decisive role. 

What is the importance of this work? First of all, we believe that the discussion 
provides strong support for the main assumptions of our model: That one can 
treat the entropy as an additional fluid, endowed with the inertial and dynamical 
properties generally associated with fluids. The entrainment between particles and 
the massless entropy plays a key role in this approach. The connection between 
this entropy entrainment and the thermal relaxation time provides an immediate 
interpretation, and illustrates the importance, of the main parameter of the model. 
The simple fact that the same mathematical framework can be used to model both 
heat conduction and finite temperature superfiuids (Andersson & Comer, 2009) is, 
in our view, clear evidence of the elegance and promise of the variational multi-fiuid 
approach. Moreover, since our model has its origin in a fully relativistic variational 
analysis, see Andersson & Comer (2007) for a review, the present discussion suggests 
a promising strategy for developing a causal relativistic model for heat conductivity. 
This is known to be a challenging problem where a number of issues remain to be 
resolved. 

NA acknowledges support from STFC via grant number PP/E001025/1. GLC acknowl- 
edges partial support from NSF via grant number PHYS-0855558. 
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